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Abstract. We study the renormalized real sample covariance matrix H = X T X/VMN- 
^/WjN with N/M — » as TV, M — > oo in this paper. And we always assume M = M(TV). 
Here X = [XjkjMxN is an M x TV real random matrix with i.i.d entries, and we as- 
sume E|TVii| 5+4 < oo with some small positive 8. The Stieltjes transform mn{z) = 
N~ 1 Tr(H — z) -1 and the linear eigenvalue statistics of H are considered. We mainly 
focus on the asymptotic expansion of E{mjv(z)} in this paper. Then for some fine test 
function, a central limit theorem for the linear eigenvalue statistics of H is established. 
We show that the variance of the limiting normal distribution coincides with the case of 
a real Wigner matrix with Gaussian entries. 

1. Introduction 

As an important branch of the Random Matrix Theory, the study towards sample 
covariance matrix traces back to the work of Hsu [8] and Wishart [IT]- The modern 
formulation of a large dimensional sample covariance matrix always indicates the matrix 
in the form of Sjy = jjX t X, where X = [Xjk]MxN is an M x TV real random matrix 
with mean zero and variance a 2 i.i.d entries and N/M — > y G (0, oo) as N, M := M(TV) 
tend to infinity. So usually the sample size M and parameter number TV are assumed to 
be with the same order. For a symmetric matrix with eigenvalues Ai, • • • , Xn, we denote 
its empirical spectral distribution by -F/v(x) = ^ YliLi l{\i<x}- And the limit of Fn{x) as 
V — t- oo is often referred to as limiting spectral distribution. In the case of Sn, it is well 
known that Marcenko and Pastur [13] firstly found: as V tends to infinity Fn(x) almost 
surely converges to the so called Marcenko-Pastur (MP) law with the density function 

1 



Py ^ = 2-nxya 2 ^( b ~ x ^ x ~ Q ) 1 {a<x<b} ; 

and there is a point mass 1 — 1/y at the origin if y > 1, where a = <r 2 (l — ^/y) 2 and 
b = a 2 (l + ^y) 2 . 

While in modern statistics, the case of N/M — > as both V and M tend to infinity 
is also very common, see El Karoui [7] for example. In this paper, we will focus on the 
theoretical aspect of such a particular case of the sample covariance matrix. It is easy to 
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see when N/M — > y = 0, the interval [a, b] will shrink to the point a 2 , so every eigenvalue 
of S will tend to a 2 then. If we centralize S by subtracting a 2 In and multiplying it with 
w M/N, the range of the eigenvalues will be enlarged to be order 1 typically. In fact, 
under the assumption N/M — > as N — > oo, it was understood long time ago that the 
matrix 

H N := A H (—X T X - a 2 I N ) (1.1) 
V N \M J K ' 

behaves similar with a Wigner matrix of dimension N on many spectral properties. 

A real Wigner matrix can be defined as a real symmetric random matrix Wjv with 
mean zero, finite variance i.i.d diagonal entries and mean zero, variance u 2 /N i.i.d above- 
diagonal entries, and all entries on or above diagonal are independent. As a cornerstone 
in the Random Matrix Theory, the so called semicircle law as the limiting spectral distri- 
bution of Wjv was firstly raised by Wigner [16] , with the density given by 

Psc ^ = V^w 2 - x?l{\ x \<2u}- (!- 2 ) 

It is not difficult to see, when we fix N and only let M tends to infinity, Hn will tend 
to a Wigner matrix with Gaussian entries under the effect of the central limit theorem. 
Thus if we let both iV and M tend to infinity with N/M — > 0, it is natural to ask whether 
Fjsr(x) of Hn tend to the semicircle law as well. A rigorous proof was given by Bai and Yin 
[5] through a moment method: under the assumption of the existence of 4-th moment of 
Xjk, Fn(x) converges almost surely to the semicircle law (|1.2p with ui = a 2 . So a question 
naturally arises, how close are the spectral properties of the Wigner matrix Wjy and the 
renormalized sample covariance matrix Hn when we set u = a 2 ? However, such a question 
makes no sense if we do not specify some particular spectrum statistics to compare. So in 
this paper, we will study some representative objects for H]y, like Stieltjes transform of 
Fn(x) and linear eigenvalue statistics of Hj\[. Then compare them to the corresponding 
results known for the Wigner matrix Wn- 

The Stieltjes transform of F^(x) is defined by 

m N (z) = [ 1 dF N (x) 



for any z = E + ir] with E G R and r/ > 0. We just consider fixed z , and specify -F/v(x) 
as the empirical spectral distribution of throughout the paper. If we introduce the 
Green function Gn{z) = (H]\r — z) , we also have 



1 1 N 

m N (z) = —TrG N {z) = j^^Gjj. 



3=1 

Here we denote Gjk as the (j,k) entry of Gn(z). Both as a powerful tool and as a rel- 
evant spectral statistic, the Stieltjes transform of the empirical spectral distribution has 
shown to be particularly important in the Random Matrix Theory. As is well known, 
the convergence of probability measure sequence is equivalent to the convergence of its 
Stieltjes transform sequence towards the corresponding transform of the limiting measure 
(for example, see [4]). So under the assumption EjAj^} < oo, one has when iV — > oo, 
rriN(z) almost surely converges to the Stieltjes transform f(z) of the semicircle law given 
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by 



-Z + Vz2_ 4ct 4 

/(*) = ^ • (1-3) 



Here the square root is specified to be the one with positive imaginary part. What's more, 
due to the basic fact |mjv(^)| < J7 -1 ) the a.s convergence of ttin(z) also implies 

f N (z):=E{m N (z)} = f(z) + o(l). (1.4) 

In this paper, we will take a step further to calculate the leading order term of the 
remainder. A corresponding work for the Wigner matrices has been taken in [IU] by Kho- 
runzhy et al., and extended by some subsequent articles to some other matrix ensembles, 
seep], [11]. Such a topic is often referred to as the asymptotic expansion of /n(z) in the 
literature. For ease of the expression and calculation, we just consider the normalized case 
of a 2 = 1 hereafter, so the semicircle law in (|1.2p will be a standard one, the general case 
is just analogous. We denote the a-th moment of X\\ by uj a and a-th cumulant by n a , 
particularly there are K2 = o~ 2 = 1 and K4 = W4 — 3 in our setting. Our first main theorem 
is on the asymptotic expansion of /n(z): 

Theorem 1.1 Consider the matrix model Hn defined in \1.1\) with a 2 = 1. And X = 
[Xjk\MxN is an M x N real random matrix with mean zero and variance 1 i.i.d entries. 
We assume M = M(N) and N/M -> as N -> 00. IfE\X u \ 5+S < 00 with some small 
positive 5 , then we have the following asymptotic expansion holds for any fixed z with 
r/>0 

f N (z) = 



+ 



Remark 1.1: By comparing with the corresponding result for the Wigner matrix in [10 |, 

we can find a different first order remainder with order rather than 1/N appears 

when M = o(N 3 ). What's more, we do only require rj > in our result, rather than 
77 > 2oj for Wigner matrix case in |10j . 

(j 1 . 3 j) and (|1.5p show the convergence of itin(z) a.s. and in expectation respectively, as 
usual, we shall consider the fluctuation of m^{z) as a further job. In this paper, we will 
deal with the more general linear eigenvalue statistic 

N 




i=i 

of Hn with the test function ip satisfies some smooth conditions. Note that iVmjv(^) is a 
linear eigenvalue statistic with (f(x) = l/(x — z). The topic of CLT for linear eigenvalue 
statistics is really classical and attractive in the Random Matrix Theory, there are a vast of 
related articles for different matrix ensembles, for examples, see [2], [3], [S], [T2J, [14] , [15] . 
However, for the large dimensional sample covariance matrices, there is no corresponding 
result on the case of N/M — > 0. So as a complement, we will discuss it in this paper. 
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For convenience, we follow the notations in [T2] to denote £° = £ — for any random 

variable £ below. And introduce the norm 

Ml = f {l + 2\k\) 2s \m\ 2 dk 1 m = e ik ^(x)dx. 

Then we can state our second result as follows: 

Theorem 1.2 Consider the matrix model defined in hl.l)) with a 2 = 1. And X = 
[Xjk\MxN is an M x N real random matrix with mean zero and variance 1 i.i.d entries. 
We assume M = M{N) and N/M ->• as N ->• oo. IfE\X n \ 5 < oo, then for any real 
test function <p satisfying ||</?||3/2+e < 00 with any e > 0, we have C^[np] converges weakly 
to the centered Gaussian distribution with variance 



V[<P] = ^ ^¥ V 2 ' S, ^ X ld X 2 




(1.5) 

Remark 1.2: Again we can compare it to the corresponding result of the Wigner matrices 
in |14| . the variance coincides with the counterpart for a Gaussian Wigner matrix (Not 
GOE!) with variance diagonal entries. This can be explained as a result of the central 
limit theorem and asymptotic independence effect for the matrix entries. We require ip 
to be real in the above Theorem, so when we deal with mjv(-z), we need to work with 
Rem ^ (z) and ImrriN(z) as in [TO] . 

Our article is organized as follows. We will provide some basic facts and tools in Section 
2, and as a warm-up, we use them to revisit the semicircle law for Gaussian case (i.e. Xu is 
Gaussian) in an average sense. In Section 3 we present the proof of Theorem 1.1, which is 
based on the preliminaries provided in Section 2 together with some extra lemmas, whose 
proofs are postponed to section 4. In Section 5, we present the proof of Theorem 1.2 in 
short and leave all technical details in the Appendix. Throughout the paper, we will use 
C and Ci(i = 1,2,3,4) to denote positive constants which may be different from line to 
line, sometimes may depend on rj. 



2. Preliminaries And Gaussian Case 

In this section we present some basic tools and facts required in the sequel, and use 
them to revisit the semicircle law proved in [S] in an average sense. For convenience, when 
there is no confusion, we will get rid of the subscript N in the notations of matrices and 
z as a variable in the notation Gn(z) below. If we set Y = \Yjk\MxN '■= {MN)~ l l 4 X, we 
have 

rp [M 

To do the job we state the following two lemmas without detailed but routine proofs. 



Lemma 2.1 (Generalized Stein's equation) 



For any real-valued random variable £ with E{|£| p+2 } < oo and complex valued function 
g(t) with continuous and bounded p + 1 derivatives, we have 

E{^(e)} = E^ i IE{5 (a) (0} + e, (2.1) 

where n a is the a-th cumulant of!;, and 

\e\ <Csup|^ +1 )(t)|E{|^r 2 }, 

t 

where the positive constant C depends on p. 

Remark 2.1. The proof of Lemma 2.1 can be found in a lot of references, for example, 
see [TO] for details. 



When £ is centered Gaussian, (|2.ip reduces to the famous Stein's equation: 

E{£ 5 (0}=E£ 2 -E{</(0}- (2.2) 

Our second lemma is on the derivatives Djk{-} with respect to the matrix entry Yjk, 
which will be used frequently in Section 2 and 3. 

Lemma 2.2 For any a,j 6 {1,2, • • • , M} and /?, k € {1, 2, • • • , N}, we have 

(i) : 

DjkGap = —{YG)j a G/3k — (YG)ji3G a k, 

(ii) : 

F>jk{YG) a p = 5 a jG/3k — G/3k(YGY T )j a — (YG)jp(YG) a k, 

(Hi): 

F>jk [Gkk (YGY T )jj] = 2G k k{YG) jk - 4G k k(YG) lk (YGY T ) jv 
where 5 a j in (ii) is the Kronecker delta function. 

Remark 2.2. The proof of Lemma 2.2 is based on the following resolvent identity for real 
symmetric matrix A and B 

(A + B- z)- 1 -(A- z)- 1 = -(A + B- zY l B(A - z)' 1 , 

which implies 

±( A + eB- zY\ =0 = -(A - z)- l B(A - z)-\ (2.3) 

Now we set A = H and A + eB = (Y + eE(j, k)) T (Y + eE(j, k)). Here E(j, k) represents 
the M x N matrix with (j, k)-th entry to be 1 and others 0. Then by \2. 3\) it is not difficult 
to get (i) of Lemma 2.2. And (ii) and (Hi) can be proved by the chain rule. We omit the 
details here, in fact it is quite similar with the counterparts in |10| and [12] . 

As a warm-up for the main task in Section 3, we use the above two lemmas to prove 
fw(z) tends to f(z) (as N — > oo) for the Gaussian case below. By the basic relation 
between a matrix and its Green function 



G = -z~ l + z- 1 GH, 
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we have 

f N (z) = -z~ l + z~ 1 ^E{TrGH} 



-z- 1 - z ^Mf N ( z ) + z -i± Y,E{Y jk (YG) jk }. (2-4) 



j,k 

Then we can use the Stein's equation (|2.2p to (|2.4p . which yields 

fN(z) = -Z' 1 - z ~ 1 \jj^fN{z) 

+ aT 1 !— = £E{G fcfc - ((YG) jk f - Gkk{YGY T )jj} 



-1 -i 



-z — z 



E{TrYG 2 Y T } - z' 1 —==E{m N {z)TrYGY T } 



-z 



-l 



4 iV + h , , _ u . [N 



yim fN(z) ~ z~\l + \J—z)E{m 2 N (z)} 



N 2K V M 
where we have used the fact that 



^4(1 + A ^zMTrG 2 }, (2.5) 



TrYGY T = Tr (l N + (\/^ + z)G\ . 

TrYG 2 Y T =Tr(G+ + z)G 2 

Note the trivial bound ||G|| < 77 for the matrix norm of G, (|2.5p implies 



f N (z) = -2" 1 - Z - 1 E{m5,( Z )} + o(Jj V 



"I 
M / 



To estimate E{m^(z)}, we need to derive a bound for Var{mAr(z)}. By using the 
Poincare inequality for Gaussian matrix entries as in Proposition 2.4 of [12], it is not 
difficult to get 

Y&x{m N {z)} =0\ ' 



AT2 



Consequently one has 



f N (z) = -z- 1 - z-'fUz) +o(^W^\ (2.6) 

On the other hand, it is well known that f(z) satisfies the equation 

/(z) = _ z -i_ z -i / 2 (z) . (2 . 7) 

By a routine comparison issue on (|2.6p and (|2.7p we have /n(z) — > f(z) as N — > oo (see 
Section 2 of [3] for example), which implies the convergence of the expected empirical 
spectral distribution to the semicircle law. 

However, for Theorem 1.1, we need to take a step further to calculate the leading order 
term of the remainder Jn( z ) — f( z ) precisely for the general distribution case. So at 
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first we need use (|2.ip instead of (|2.2p . As a result, more involved estimate towards the 
derivatives is required. What's more, the Poincare inequality is no longer valid, we need 
to estimate Var{mAr(z)} by a martingale difference method as in |14| . 



3. Asymptotic Expansion For /n(z) 

To prove Theorem 1.1, we begin with the basic idea raised in [10], use (|2.ip to (j2.4j) 
and estimate the derivatives. However, the standard process will be too involved for our 
matrix model . As is shown in [lOj , for the Wigner matrix, every term in the expansion 
formula can be factorized into the entries of Gn{z) which can be bounded by 77 -1 . Owing 
to such a trivial bound, many estimates hold obviously, especially for the remainder term. 
And the approach of iterating the expansion process also works well. However as we will 
see, for our case, the derivatives are in more complicated forms. There are no trivial bound 
for some factors of the terms (see (|3,8|) . (|3.9|) 1. as well, iterating the expansion process will 
bring factors with new types. And to bound the remainder, we need to use a truncation 
technic at first. 

Now we truncate Xj k at r := (MA) 1//4_ * with some small positive t < (say). Then 
we denote the empirical spectral distribution of the matrix defined in (|1.1|) with truncated 
entries Xj k := Xjkl^ Xjk \< T ] by F^{x). Further we denote the stieltjes transform of F]^(x) 
by m T N {z) and set f^{z) = Em T N {z). By taking into account that \rriN{z)\, \m T N (z)\ < T]^ 1 
and 

^2^{\X jk \ >t}< CMN-t-^ +&) < C(MN)~^~ t 

jk 

we have 

\f N ( Z ) - r N ( Z )\ , 

After the truncation, all of the first 5 moments of Xj k will be modified. For example, 
\E{X jk l {[x { < T} }\ = |E{X ife l { |x, fc |>r}}| < r-^E\X jk \ 5+s < C{MN)- 





-l-t 



Similarly we can calculate the modification of the a-th moments of Xj k under the 
truncation for all a = 2,3, 4, 5. So without loss of generality, we can always assume 
Xj,k, (j, k) e {!,-■■ , M} x {1, • • • , N} are i.i.d and 



\X jk \ < (MA0 1/4 ~*, \E{X jk }\ < CiMN)- 1 -*, 
E{A|J = 1 + C(MN)~i-\ E{X« k } =u a + o(l), a = 3, 4, 5 
(jk\ 



E{\X jk \ 5+s } < C. (3.1) 



And for simplicity of the notations, we still use Jn(z) to denote fj^{z) below. What's 
more, we need to present two lemmas before the rigorous proof of Theorem 1.1. The first 
lemma is on the estimate of Var{mAr(z)}, as we have seen in Section 2, it is a necessary 
ingredient for the final result. The second lemma is on an estimate towards every diagonal 
entry of the Green function, which can help to overcome effectively the difficulties we have 
mentioned above. Their proofs will be postponed to Section 4. 
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Lemma 3.1: Under the assumptions of uj\ = and U4 < 00, we have 



for any fixed z with r] > 0. 



V&i{m N (z)} = 0[ 



To present the second lemma, we need to introduce the matrix 

-1 



G = [Gjk]MxM 



YY 



T 




Lemma 3.2: Under the assumptions of Lemma 3.1, for any fixed z with r/ > we have 
the following two estimates 
(i): for any k G {1, 2, • • ■ , N} 

2 



E 



Gkk + 



1 



(ii): for any j G {1, 2, • • • , M} 



E 



G 



j.i 



M+z + f N (z) 



<r( N )3 + r N 

<<?!(-) +C 2W 



(3.2) 



(3.3) 



Remark 3.1 In fact, we can replace (z + /^(z)) -1 in (|3.2p by /jv (z) or even f(z) with the 
same bound on the right hand side, but as a lemma, we shall not do that before we get 
the final conclusion on the differences between {z + /^(-z)) -1 , /at(-z) and f(z). 



Proof of Theorem 1.1. Use (|2,ip to (|2.4p . we have the following expansion for the general 
distribution case 

f f \ -1-1 ft*, f V 

+z~^ £ t^it E ^r^(yc), fc + ^, (3.4) 



^ (MiV) ; 



where represents the a-th derivative with respect to Y a p and 

Rn < C±—^-£supE jk \D%(YG) jk \. 

Observe that in (13.41). the o = term 



Ut7E{x u } J]E(y% = °(v ^ v ^) 



iV(MiV)V4 

by taking into account (|3.1|) and the estimation (|3.13|) proved below. So we will focus on 
the estimations of a > 1 terms and i?/v in the sequel. Here sup jfc means the supremum 
is taken w.r.t Yj k and Ej k represents the conditional expectation w.r.t Yj k . By using 
Lemma 2.2 repeatedly we have 



D jk (YG) jk = G kk -[(YG) jk } 2 -G kk {YGY T \ 



33 
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D 2 {YG) ]k = -QG kk {YG) jk + QG kk {YG) jk (YGY T ) lj +2[{YG) jk f (3.5) 



D%(YG) jk = -Q{G kk ) 2 + mG kk [{YG) ]k \ 2 

+12(G kk ) 2 (YGY T ) jj - ZQG kk [{YG) ]k ] 2 {YGY T ) j] 
-6(G kk ) 2 [(YGY T ) n ] 2 - 6[(YG)^ (3.6) 
D%{YG) jk = l20[G kk ] 2 (YG) :jk - 240G kk [(YG), k } 3 

-24Q[G kk ] 2 (YG) jk (YGY T ) jj + 2A0G kk [(YG) lk ] 3 (YGY T ) lj 
+120[G kk ] 2 (YG) jk [(YGY T ) jj } 2 + 2A[(YG) jk f 

for any (j, k) G {1, • • • , M} x {1, • • • , N}. Observe that for any integer a > 1, Dj k (YG)j k 
is a linear combination of terms in the form of 

^(01,02,03) := [G kk r[(YG) 3k r[(YGY T ) jj r (3.7) 

with nonnegative integers ai, 02, 03 satisfying 01+02+03 < a+1 . Moreover, by Lemma 2.2, 
it is easy to see one (YGY T )jj must appears together with one G kk , thus we always have 
a\ > 03. By observing that 



T ~ T , /M , ~ 

ycy T = cyy T = i M + (W — + *)g 



and the trivial fact ||G|| < r/ 1 we have 



<M 
~N 

Similarly we have the following estimate 




|(ycy T y < c x /^. (3.8) 



/ \ 1/2 /m\ 1 / 4 

\(YG) jk \< r£\(YG) jk \ 2 j =\(YGG*Y T ) JJ \ 1 / 2 <C^-j . (3.9) 

But these two bounds are too bad, thus we need to figure out the higher order term Rjy 
more carefully. Fortunately, when we take expectations in (|3.4|) . we can use the following 
estimate 



E 



(YGY T )jj 



In(z) 



% + z + f N (z) 



^ c ^ )2 + c *w (3 - 10) 



which is a direct consequence of (|3.3p and ([37 
To deal with Rn, we observe that 

1*^(01,02,03)1 < C\(YG) jk n(YGY 



T\ . 103 
)33\ 



C| E YjaG ak \ a2 1 YjpYj^G 



/3-r I 



a a ' ^ /3, 7 ,9,7 

< CiV^Y/j^. (3.11) 
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In the second inequality we used Cauchy Schwartz inequality and in the last step we used 
the fact that 

\G*j\ 2 = [G*G] j3 <C, Yl \°^\ 2 = TrG * G = CN 



Pa 



E Y h Y i 

Pa 

So by (|3.1ip we have 



2 

31 



E^ 



a 2 +2a 3 

supE ifc |F jfc (oi,oa,o 3 )| < CN% ( sup \Y jk \ a2+2a * + E( Yf a ) * ) 



a 2 + 2a 3 
4 



Then it is not difficult to see 



\R N \ = o(\ — V -4" 
11 V M AT' 



by taking into account d2 + 2a3 < 5. 

Now we can turn our attention to (|3.5 j) and (|3.6p . Firstly, we shall provide some coarse 
estimate on Yljk \0^G)jk\ a2 ■ When 02 > 2 we have 



^[(G*y T ) fc ,(yc) 



02 
2 



j,A: 



E 



G fcfc + (y^ + z)(G*G) fcfe 



£12. 

2 



= 0(M*iV » ). 
When 02 = 1, we can use the elementary inequality 



(3.12) 
(3.13) 



j,k V jfc 

And for this time we only need to bound the expectations, so (|3.10p can be used. By 
Lemma 3.2, (|3.10p . (|3.12p and (|3.13p . we can discard all the terms except for 

-6F jfc (l,l,0), -6^(1,0,0) 

from ()3.5p . ()3.6p respectively. For F^(l, 1, 0), according to Lemma 3.2 and (|3. 13[) . it 
suffices to estimate Yljk ^{(YG)jk} instead. To get a lower bound than the coarse one 
(|3.13p . we need to iterate the expansion process, the main term of Yljk^ , {( Y G)jk} is 
enough for us. We do it as follows: 

E^2(YG) jk = ^E(Y ja G ak ) 

j,k j,a,k 
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E 

j,a,k 



EDj a G ak + — 



where 



ejvl < C T7Tf snpE ja \D ja G ak \. 



(MiV)t 



+ e N , (3.14) 



MN 



J a 



By using Lemma 2.2 again we can get 



D jaG ak 



.1 

Dj a G a k 



— (YG)j a G ak — (YG)j k G aa 
—2G aa G a k + 2G ak [(YG)j a ] 2 
+2G aa G a k(YGY T )jj + 4:G aa (YG)j a (YG)jk 



(3.15) 
(3.16) 



ak — 18G aa G a k(YG)j a + 6G 2 ta (YG)j k 

— l8G aa G a k(YG)j a (YGY T )jj — 6G ak [(YG)j a ] 3 



-18G aa (YG) jk [(YG) ja ] 2 - 6[G aa ] 2 (YG) ]k (YGY T ] 
Inserting (|3.15[) and (|3.16p into (|3.14|) one finds that 



33- 



3,k 



I T AT ^-^ 



y/MN 



3,k 



+ - 



(YG 2 ) jk + (YG) jk TrG 

MG aa G ak + G ak {GY T YG) c 



G aa G ak Tr(YGY T ) + 2{GY T YG) ak G, 
Similarly to (|3.13p we can provide 



EY J (YG 2 ) jk = 0( K (MNY< 

j,k 

The other terms in the expansion are easier to estimate by a similar calculation as we have 
done above. Moreover, similar to the case of Rn it is not difficult to check that e^r has no 
contribution to the main term of ^ - k E{(YG)j k }. So we can get a bound as 

EY / (YG) Jk = o^nm 

Therefore, under the estimates above, we arrive at 



f N (z) 



-z — z 



-1! 

N 



l ^Y,[(YG) jk f 

3,k 

EY,G kk {YGY T ) 3J - z^^Ej^Gl, + o{±). 

j,k k 



By ()2.5p and Lemma 3.1, we can rewrite it as 



!n{z) 



-z ^—E{TrG 2 (z)} 



N 2 
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which implies 



In(z) + 




(3.17) 



Z + fN(z) 

Furthermore, we know both Jn(z) and f(z) are analytic function of z. Then by (jl.4p 
and the Cauchy's formula, we also have 



d 

dz 



(/*(*)) 



d_ 

dz 



(f(z)) + o(l). 



Now we note that 



^TrG\z) = = Tz if{z))+0{1 

Together with (|3.2p we finally have 

r W(f(z) + zf(z)) 



z + 2f(z] 



+ o(l). 



In(z) 



f(z) 

1 



+ 



N 



M z + 2f(z) 
f(z) n 4 f 2 (z) 



(z + 2f(z)y z + 2f(z) 



N 

+ 0(-) + o 



This concludes the proof by inserting the basic relation z + 2f(z) = f(z) — 1/ f(z). □ 

4. Variance Estimates 

In this section, we will state the proofs of Lemma 3.1 and 3.2. In fact, we can use 
the asymptotic expansion method again to estimate Ysa{m^(z)}, but it is really more 
complicated than in the Wigner case. We will use a martingale difference method used in 
|14| very recently. For convenience we introduce the notation 77V to replace TrG below. 



Proof of Lemma 3.1. If we denote E<&{-} and Efc{-} as the expectation w.r.t the ran- 
dom variables of the first k columns and A;-th column of Ijy respectively, by the classical 
martingale method in [6J, one has 

N 

Var{ lN } = Y, E {\ E <k-iilN}-^<k{lN}\ 2 } 
k=l 

N 



^E{|E< fc _ 1 { 7 j V -E fc { 77V }}| 2 } 



k=l 
N 



< 



^E{| 77 v-E fc { 77V }| 2 }. 



(4.1) 



k=l 



Now we let y k be the k-th column of Y/v and B k the M x (N — 1) matrix consisting of the 
other N — 1 columns of Yjy, then we have 



H 



yi ■ yi 



Bfyi 
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With the notations 



M 



= BlBk — \ —In~i 



N' 



and = {H^ - z)' 1 , we have 

l + yi -B x Bj ( Bl BT-J§. 
TrG - TrG (1) = ^ 



z yi 



//. ■ Hi. - \l t- - * ~ Vi ■ B X B{ ( B x Bj -JM-z\ yi 



l + U 
V 



(4.2) 



and Gn = V . Here we used the basic identity 

Bi (b t iBi - ^ - *) ~ > = B lB f (B X B? -J§~z 
To estimate ()4.ip . we only need to deal with the first term 



n = 1,2. 



E{| 7J v-Ei{ 77 v}| 2 } =E 



1 + Z7 



V 



Ei 



l + U 
V 



(4.3) 



since the others are analogous. So it suffices to estimate E^UV' 1 - EijW -1 }! 2 } and 
E{|1/ _1 — Ei{l/ _1 }| 2 }. We will only present the estimate for the first one below. Clearly 
we have 



Ei 



U-E X {U} V-E^VjU 



Ei{V} 



Ei{y} v 



By observing that |Ei{V}| 1 < rj 1 and 

y 1 . Bl Bf\B 1 Bf-^- 



\U/V\ < \U/lmV\ 



A yi 



< 



n + rjyi- B x Bl\B x B\ -J§- z\~ 2 yi 



it suffices to provide the following estimates 



Ei{|C/-Ei{C/}| 2 }, Ei{|y-Ei{V}| 2 } = 



N 



For simplicity we introduce 



B, BT I BlBf — A / -r-r — z 



' M 



k \MxM '■= n\r>\ D\D X 



N 



(4.4) 



(4.5) 



'M 



Ml 2 ' = [M% lxM := Bi-Bf ( B lB f -\ljj-z 



Thus we have 



U-E 1 {U} = Y^Mlf^i + Y,^^-^^' 

i^j i 
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V-MV} = yx-yt-^-^M^YaY^-Y^M^iYl-E^Yl}). 

ijdj i 

It is not difficult to get 

E 1 {\U-E 1 {U}\ 2 } < Cj ^^\ M lf\ 2 = C J^ Tr \ M[2] \ 2 ' ( 46 ) 

E X {\V - E 1 {V}\ 2 } < ClJ ^Tr\M^\ 2 + C 2 i. (4.7) 

Now we denote eigenvalues of by fj,^ ■ ■ ■ fJ^]_i- Note that ^a\a = 1, • • • , N — 1) are 
also the eigenvalues of the matrix 

H {k) := B k Bl - S-^Im, 

which has a (M — N + l)-multiple eigenvalue —\J~^- So it follows 

N-l (J- 1 ) _L fM)2 

Tr|MN| 2 = £ ^ V"' = „ = lj2> 

l/i* * *y I 2 71 

a=l I pa z l 

which implies (|4.5p by taking into account (|4.6p and (|4.7p . □ 

Proof of Lemma 3.2. We will only prove for G\\ and Gn, the others are analogous. As 
we have shown above 

Gn = V- 1 = =i . (4.8) 

VI- Vi- ^-z-m-MKyx 

We denote the unit eigenvector of corresponding to the eigenvalue fj,^ by Va = 
(v { a\l),--- ,v { a\M)) (a = l,--- ,N- 1) and set 

It is easy to check Ei{^ 1} } = 1. Then by (1431) . we have 

a u = ! 

/ M _ _ l v-^ -1 _ V " iv ■ ,Sa 
Vl ■ Vl y N z /aJjV ^a=l ~pT 



Ha ~Z 



-z - mjv(z) + ri 



where 



iV-l + 1 

iV ^ ' 

a=l 



/ / AT— 1 / Jv , 

(i/i • yi - a/ -T7) + (m N (z) - - ^ 



a =l Po ^ 



"iv S -tit— 

a=l P« Z 
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m N (z) 



1 V A/ "I" 1 



iV 



E 



1 pa — 2 



+ y-Ei{F}. 



As a consequence we have 



1 



^ + In(z) 



f N (z) - m N (z) + n 



(-z - mAr(z) + ri)(z + 
< C\{f N (z)-m N {z))+ ri \. 



Thus we obtain 



E 



Gn + 



< C(Var{m N {z)} +E 



miv 



,yW z 

=1 Ma — z 



+ E|V-Ei{V}| : 



Note the trivial bound 




1 N - 1 



N ^ _ z 

1 /-*a 2 




(4.9) 



(4.10) 



and the fact that 



m N (z) 



1 1 

iV ^ (i) _ 

a=l Z 



dF N {x) 



x — z 



x — z 



1 

N 
1 

< — 



NF N (x)-(N -l)F^(x) 



{x — z) 2 



dx 



dx 



7T 



(4.11) 



j (x- z y 

where (x) is the empirical spectral distribution of ijW . In the inequality above, we 
use the well known interlacing property between eigenvalues of an Hermitian matrix and 
its submatrix as: 

Ai<^S 1) <A 2 <---<^ ) _ 1 <A Jv . 

Combining Lemma 3.1, (|4.5p . (|4.10p and (|4.1ip we conclude the proof of (i) of Lemma 3.2. 
It is similar to get (|3.3p for Gn, we now turn to 



H = YY 



r 




M- 



If we denote the /c-th row of Yjv by and the (M — 1) X N matrix consisting of the other 
M — 1 rows of Y/v by D k . And further introduce the matrix 



H^ = [H%% xN =:D T k D k - x r-I N 



' M 
~N' 
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(k) 

with its eigenvalues fl a , a = 1, • • • ,N. Then we have the following representation 



H 



yi ■ yi - y w yi^i 

D 1 yJ D X D\- Mlu-i 



If we use the E\{-} to denote the expectation with respect to yi, and set 

\ -l 



We have 



V = mDjD 1 D{Di - 




— z 



Vi ■ 



G 



ii 



V 



~v w ~ z ~ m N{z) + fx 



where 



h = yi-yi+m N (z) + V 

= y 1 .y 1 + mN (z)--J2 1 ^m + V-E 1 {V}. 



a=l H'Ct Z 



Now we denote the set of event 



=: { \/n(z) - m N (z) + fi(z)\ > 




and observe that on Cl^ there exists 




' '- + »'..V( -) - '!(-))( \ljjr + Z + f N (z 



Then we have 



E 



E 



Gn + 



f + z + f N (z) 



f N (z) - m N (z) + fi(z) 



C 



( y/f + z + mjv(z ) - fl (z))( y/f + z + 

El/jvOO-mjvO^+riO*)! 2 , 



<f> 2 + >°> 



where we used the fact that 




+ z + m N (z) - fi(z] 




+ z + f N (z) 



> C 




(4.12) 



(4.13) 



holds on the full set ft. To see this, We denote the unit eigenvector of H^> corresponding 
to the eigenvalue by Va^ = (ia^ (1), ■ ■ ■ , Va^ {N)) (a = 1, ■ ■ ■ ,N) and set 

f£ 1 ) = v / MiV|yi-^ 1 )| 2 . 
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Then we have 




+ z + m N (z) - f\{z) 




+ z-yi 



N 
M 

— + z+ , 

N y/MN 



+ 




Vi 




N 

£ 

a=X 



M , ~(1) 
jY \ fa 



i e 



(i) 



N 



a=l fo 



?(1) 



Taking into account that /} 



(i) 



> 



ReS > -C\ —TmS, 



+ «)(1 + S). 
-y and fi 1 ^ > 0, we have 



Im5 > 0, 



which implies 



N 



Then (|4.13p is a direct consequence of (|4.14p and (|4.15p . 
Now we proceed to the estimate of (|4,12p . By observing 



P(0 ) < ^E\f N (z)-m N (z)+h(z)\ 2 , 



it suffices to provide 



N 



1 



nf N (z) - m N (z) + h(z)\ 2 < Ci— + (7 2 - ; 
which can be derived similarly as what we have done to ()4.9p . 



(4.14) 



(4.15) 



□ 



5. CLT For Linear Eigenvalue Statistics 

To prove Theorem 1.2, we will follow the recent article [2] by Shcherbina. For there 
are only some technical differences, we will only state the main body of the proof in this 
section, and left all the technical details to the Appendix. 

Proof of Theorem 1.2. Similar to [T3], as we will see, cog is needed in our proof. So firstly 
we need to truncate the random variable Xjk at (MAT) 1 / 4 , and then re-centralize it. To 
use the truncated matrix, it is necessary to show its linear eigenvalue statistics have the 
same limit distribution as the original one at first. If we denote the truncated matrix as 



Y 
Y 



[Yij]MxN, 
[Yij]MxN, 



y. 

1 7 7 



^i^-fliijl^l}' 



Y — Y 



And further introduce 



YY -h ,N ' 



C N [(f\ = Tnp(H), 
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H N = Y T Y - J ^I N , C N [<p] = Tr<p(H), 

H N (t) = [Y + t(Y - Y)f[Y + t(Y - ?)] - ^-I N . 

We denote the eigenvalues and corresponding eigenvectors of H^(t) by Aj(t) and Vi(t) 
(i = 1, • ■ ■ , N) below. It suffices to prove that 

E {| e ia^M _ < 2F{Y ^Y} + \x\\E{£ N [cp]} - E{£ N [tp]}\ -> 0, (5.1) 

E || e i^M _ e «^M|} < 2 |x|E{|L i vM - £n[<p]\] -»• 0. (5.2) 

as ./V tends to infinity. The proof of (|5.2|) is similar to (|5.ip , we only give the proof of (|5.ip 
below. P{y ^ y} is obviously o(l) by the truncation. Also, we have 

E{\C N [ip]-C N [<p}\} 

= f dtE{j2\v'(\(t))m\} 

Jo 

< Halloo f 1 dm^vfm'itMt)} 

Jo 

<y\\oc [ dtE{Tr\(Y -Y) T Y + Y T (Y -Y) + 2t(Y -Y) T (Y -Y)\} 
Jo 

< Halloo f dtE{J2 \[(Y - YfY + Y T (Y -Y) + 2t(Y - Y) T (Y - Y)} %] \} 

< 2|| ¥ >'|| 00 f dty2(E{\(Y ki - Y ki )\ ■ \Y kj \} + iE{|(y fc i - Y ki )(Y k j - Y k j)\}) 

< 2\\ (f /\\ o0 [ dt{ Y\E{\Y kl - Y ki \ ■ \Y kl \} + tE{\Y kl - Y ki \ 2 }} 

Jo 1 ~k 

+ y, YJ^\ Y ki - ■ E ftkj\ + m \ Y ^ - ?ki\ ■ w kj - ? kj \}\. 



(5.3) 



i^j k 

By the truncation, we have 

E\Y kj \ < C(MN)-\, 

E\Y kl - Y ki \ < (MN)-\ J \X ki \l {lXktlHMN)i} dP < C(MN)-l. 
E\Y kl -Y k ^<{MN)-\ J \X kt \H {iXk ^ {MN)h dP<C(MN)-l 

Then it is easy to check (|5.3j) tends to as A'" tends to infinity, so (|5.ip holds. 
So without loss of generality, we may and do assume in the sequel 

E{X jk } = 0, E{X 2 jk } = l + o(l), E{Xf k }=co 4 + o(l), (5.4) 

E{X^ k } < C{MN)*, E{X s jk } < C(MN)i. (5.5) 
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Relying on (|5.4j) and (|5.5p we have the following lemma which collects all estimates we 
need to derive the central limit theorem for £jy[y>]. 

Lemma 5.1 Under the assumptions |5.^[ ), 15. 5\) and owr basic assumption N/M — > 
as N — > oo, we have for any fixed z with rj > and Z\,Z2 '■ Imzi^ > 0, 1 > 5 > the 
following estimates hold: 

N 

Var{ 7w } < CN- 1 ^{(G^z)^})^- 3 - 5 + ?T 3+5 ), (5.6) 

i=l 

E{(V°) 3 }, E{{U°) 3 }, E{|y°| 4 }, E{\U°\ A } = 0(N-'i), (5.7) 
NE 1 {V°(z 1 )V°(z 2 )} = w 4 - 1 - ^^Tr^^) + MW(z 2 )] 

+^(V / f + Z ^\fw + ^ {1) (ft)7^ (1) (^), (5-8) 
VarfiVE^F ^)^ ^)}}, Var{iYE 1 {y°(z 1 )C/°(z 2 )}} =o(%-V±=), (5.9) 



IE{|tXt| 4 } < CN l / 2 r 12 + o(iV 1/2 ), E{| 7 ^ (1) - 7 ^| 4 } = 0(iV- 3 / 2 ), (5.10) 

- f(z)\, \E- l {V} + f(z)\ = o(l). (5.11) 

We postpone the proof of Lemma 5.1 to the Appendix. Here we need to mention that 
the constants in (|5.6p or (|5.10p are independent with rj. 

By Proposition 1, Proposition 3 of |14j and (|5.6p we only need to prove Theorem 2 
for the functions (p = (p„ which are the convolution of some </?o with the Poisson kernel 
P v = n/i:{x 2 + n 2 ). And tpo is restricted to satisfy f \tpo(\)\d\ < C with some positive 
constant C. For such test function <p we have 

£n[<p] = ~ [ ¥'o(A i ) Im 7Ar(^) d / i ' z )1 = fi + in. 



7T _ 

Follow the notations of |14j . we set 

Z N {x) = E{e ixC »M}, e(x) = e ixC °^\ e x {x) = e ix ( C N } M)\ (5.12) 

where £^[99] stands for the corresponding linear eigenvalue statistics of Observe 
that Theorem 2 can be proved by providing that 

^-Z N (x) = -xV[<po,rj\Z N (x) + o(l). (5.13) 
ax 

To do this, we introduce 

Y N (z,x) =: E{TrG(z)e°(x)} = ^ E{G kk (z)e° (x)}. (5.14) 



Thus we have 

dx JV v ' 2ir 



d Z N (x) = - 1 - / cp (fi)(Y(z ft1 x) -Y(z^,x))dfi. 
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We only need to deal with the first term in the summation (|5.14p . the others are 
analogous, so we can deal with N'E{V~ 1 e°(x)} := T\ + T 2 instead, where 

Tx = -NE{(V- 1 )°e 1 (x)}, 

T 2 = -NE{(V- 1 )°(e(x)-e l (x))}. 



For the calculations towards T\ and T 2 are similar to those in the case studied in |14j . we 
will not state the tedious process. Indeed, by inserting the estimate in Lemma 5.1, we 
easily obtain 

Tx = f z (z)Y N (z, x) + Oliv.ffl 



N V M ) ' 



and 



where 



To = ixZ 



N 



(x) J dfj,(p (fJ,) 



D(z,z^) - D(z,zJ | 



2ivr 



D(z,z.) := 2/ 2 (z)/ 2 (z M )(l + /'(^)) 



/(-*) ~ , ^4-1 



z - z,, 



d ff(z)-f(z fl ) 



z - z„ 



So we have 



/D(z,z 



D(z,z lx ) - D{z,z,j) 
P{z)) 



Thus the variance in (|5.13|) can be represented by 

1 

4^2 



V[<po,rj\ 



¥>o(Mi)¥>o(^2)( C(z lxl ,z IJa ) + C( 



where 



C(z,Zfj,) 



D(z,Za) 



(5.15) 



i-f 2 (zy 

We can get (|1.5p by comparing (|5.15|) to (2.40) of [H], then we complete the proof. □ 



6. Appendix 

Proof of Lemma 5.1. We begin with (|5.6p . we will use (|4.ip and (|4.3p . and so we need to 
estimate (|4.4p more carefully. According to (|4.6p and (|4.7p . we have 

2> 



Ei 



y-Ei{y} 



<Ci 



1 



+ C 2 



1 1 



MN \z + ^=TrM^\ 2 N \Ei{V}\ 2 

1 VMA 1 



(6.1) 



and 



E-, 



C/-Ei{C/} 



Ei{V} 



C 



rr|Mt 2 lp 



MA | z + TrMWl 2 ' 

1 VMA 1 
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For (pUj) . if we set N^ = BfBi, we have 

TrAfW = TrG^N^, Tr\M^\ 2 = TrG^N^G^*N^\ 



So 



T^ Tr \ Mll] \ 2 
MN 1 1 

= i TrG m N m G w* N m 

MN 

< _^-[TrG'( 1 )ivWG( 1 )*] 1 - <5 [rrG«(iV( 1 )) 1 ^G( 1 )*] <5 
- MN 



< — 

- iV 

1 

< — 

- iV 



1 



1 



1 r 1 / N \ 



i+s 



1 5 



n\mJ v ; 



Furthermore we have 



1-8 



N~ 



M 



-2d 



-25 



which implies 



E 



ImTrG (1) iV (1) = 7?TrG (1) iV (1) G (1) *, 
1 7V|A/W| 2 



MN 



z + 



'MN ' 



< £,-«E 



1+6 



|E!{F}| 



1+5 



1 

- N 1 



-l-Si 



1+1 

LxWxM) 6 l*Ei|G n | 1+i 



(6.2) 



Here we have used Jensen's inequality |Ei{V}| 1 < Ei{|V| 1 }. 

If we control (|6.2p by using Theorem 5.9 of [1] by replacing y with N/M, we can get 

2 ) 1 

< C— E|Gii| 1+ *(jT 1 - <5 + *T 1+ ^ 



E-, 



y-Ei{^} 



i 

iV" 



Ei{F} 

Then by the fact that ImV = —r]U — rj, we also have 



n\ v w ^ } l\ 2 } < cl(n\Gn\ 1+5 })(v- 3 - 5 +v- 3+s ) 



Ei{y} v 



and 



r ii| 



})(r? -3-5 + r? -3+ 5) 



Next we turn to the proof of the first inequality of (|5.10|) . which will be used in the 
proof of (|5.7p . We use the following inequality for martingale (see [6]) 



N 



k=l 
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Similar to (|4,3p . it suffices to check 



and 



E{\U - Ei{C/}| 4 } < CN-27]- 8 + o(N~ 



E{\V - Ei{y}| 4 } < CN~2r}- 4 + o(N~ 



We only present the proof of (|6.4f) below, (|6.3[) is analogous. Observing that 1 
(BxG^Bj^i, $£E§ and (IBTTTJj) are valid as well. We have 

,M 
~N 



which implies 

E|T/-Ei{y}| 4 



< 



+ 



w 8 



+G377^7wE{(rr|M| 2 ) 2 } + C4- 1 



7^8 



' (MN) 2 LV 1 1 ' J MiV 2 
< CiV"ir ? - 4 + o(iV"t). 
For (|5.7j) . we only deal with E{|V°| 4 } below, the others are similar. 

y° = y-Ei{F} + Ei{y}-E{y} 

= V - Ei{V> + , 1 (Tr flf M - F.Tr fl/f W) 

So use the first inequality of (|5.10j) to 7^ , together with ()6.5[) we can get (|5.7p . 
For f|5.8|) . note the following expansion: 

iVE^y ^)^ ^)} 



(W4 " 1) 



CJ4 — 1 

M 



Tr 



+ 



tJ4 ~ 3 VAfM^ Mi-fM/ 1 * 2 
M 



Y J M^{z 1 )M^{z 2 ) + -Tr[M^(z 1 )M^(z 2 )} 



~[(TrM [1] ( Zl ) -ETrM [1] ( Zl ))(TrM [1] (z 2 ) -ETrM [1] (z 2 ))] 



(CJ 4 " 1) 



UJ4 — 1 

M 



Tr[M^( Zl ) + M^(z 2 )} 



W4 — 3 



M 



M 
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To deal with the first estimate in (|5.9|) . we only need to take care of the variance of the 
third and fourth terms of (|6,5p . Using (|3.1U|) to again we can get 

Var{M|l}<C l( ^ + C 2 i 
together with the trivial bound |-M^| < C \J~~^f we have 

i 

i 

^ iE E (^ ] ^)^ l0 ^) + ^f l0 (^i)^ ] (z 2 ) 

i ^ 

< Ci— + c 2 — . 

- M N 



Also we have 



TrM [1 \z 1 )M [1 \z 2 ) 

_ ^ Off + ? 
2-*> i (i) M (i) \ 



M M Z1+Z2 , , 



-, N-l 



^ C (!) \/ (1) \ 

a=l (/4 ~ 2l)(/-4 - 



JV-1 



i Zi+ Z2-, ^1 — ^2 



a=l (fia ~ Zl){^ - Z 2 ) 



+ N- 1 



M / M Zi + Z 2 w \ 

lv " ZlZ2 + ( V iV + + Z2) 



Tr 



G (i)( Zl )_ G (i)( 22 ) 



2l - ^2 



+ty f + ^)Tr(G«(, 1 ) + G«(, 2 )) + AT - 1. 
Using ()5.6p to TrG^ = we have 

The estimate towards Y&r{NE,i{V° (zi)U° (z 2 )}}) is straightforward by taking into account 
(G«) 2 = dG^/dz. 

The second estimate of (|5.10p is a direct consequence of (|4.2p and (|5.7p . The first part of 
(|5,11|) is just the consequence of (|1,4|) if we replace Jn(z) = Jn/N by jff /N, and second 
one follows directly from the fact 



m} = z + ^7? + 0(J ^) = , + /( : ) + o( I ) 
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and 



/(*) = - 



1 



* + /(*)' 



So we complete the proof. 



□ 
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